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ABSTRACT 


Molecular study is important to detect variations and similarities among species from the same genus, in case 
if they do not encompass any morphological or physiological differences. The study was conducted to 
differentiate among species of Meretrix spp. (Meretrix lyrata, M. meretrix, and M. lusoria) obtained from two 
locations in Malaysia through the phylogenetic tree. The adductor muscle tissues were used to extract DNA and 
to perform other procedures; the samples were subjected to analyses using PCR and gel electrophoresis. The 
multiple sequence compatison was conducted by MUSCLE and the phylogenetic relationships were established 
using Maximum Likelihood (ML) statistical methods with MEGA 6.0 statistical software. M. /yrata samples 
showed 99% similarity to the three accessions sequence, where M. Arata indicated 87% similarities, and M. 
meretrix showed not more than 89% similarities from the deposited sequence. The nucleotide base composition 
sequences consisted of the mean of Thiamine (T) 37.9%, Cytosine (C) 15.4%, Adenine (A) 27.4%, and Guanine 
(G) 19.4%. Maximum Likelihood (ML) analysis was conducted using the Tamura 3-parameter model to establish 
five major clades on Meretrix spp. and two out-groups clades significantly different from the Meretrix spp. These 
major clades were closely related to each other at the 50% evidence of bootstrap, which grouped as genus 
Meretrix. The present study on Meretrix spp. from the Sarawak locality was able to differentiate COI sequences 
between M. /yrata, M. meretrix, and M. lusoria. M. lusoria was close related to M. meretrix with strong bootstrap 
supporting evidence at 96% scoring. Moreover, M. /yrata was inferred as the ancestor to M. meretrix, and M. lusoria 
from Sarawak, Malaysia. 
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INTRODUCTION 


DOE: 10.11598/btb.2022.29.3.1583 


Genetic variation can be characterized into 
two groups which are intraspecific and 
interspecific (Arruda e¢ a/ 2009; Ehlers et a/. 
2016; Layton e¢ al 2016; Thia et ai 2016). 
Intraspecific is the variation within individuals 
from similar species such as subspecies (Stevens 
et al, 2010). On the other hand, interspecific is 
the variation within individuals above the 
species level. Genetic variation of natural 
populations is thought to be governed by the 
mutual effects of irregular genetic drift, 
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restricted gene flow, and variance of selection 
pressures. Intense natural selection is an 
important factor influencing variation among, 
and within groups (Endler 1986). The natural 
selection can be stabilized the species characters 
such that genetic variation within the population 
might be reduced. Disruptive selection may lead 
to genetic divergence, and therefore, may 
eventually increase genetic variation among 
population (Backeljau e¢ a/, 2001). 

Mitochondrial DNA (mtDNA) is one of the 
DNA markers that are popular for DNA 
identification. This marker is the maternal 
inheritance that supplies information absent in 
nuclear markers (Okumut & Ciftci 2003). 
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Therefore, selecting tissue from the female 
individual is essential to ensure information 
regarding the species ancestor. However, the 
maternal inheritance by mtDNA _ can_ be 
restricted by the geographical with no mixing of 
mtDNA _ haplogroups from the different 
geological regions (Luttikhuizen e¢ a/ 2003; 
Mishmar e¢ a/ 2003; Layton ef ai 2016). 
Furthermore, inheritance characteristics on a 
particular defect gene could be recognized using 
an mtDNA marker (Wallace et ai 1999). This is 
significant, particularly in aquaculture, to ensure 
cultured species are free from any inherited 
diseases from the parents. Moreover, mtDNA 
could be applied for species identification and 
phylogeny studies due to the particular inherited 
gene variance. 

The typical mitochondria gene used for the 
phylogenetic study is cytochrome C oxidase 
subunit I (COI). This mtDNA marker has been 
applied for genetic relationship and phylogenetic 
analysis for Meretrix spp. from Japan, Korea, and 
China (Yamakawa ef a/, 2008, Chen e¢ a/ 2009; 
Torii e¢ af 2010; Kim & Yoon 2014; Sato ef al 
2016). Parentage assignment from a_ similar 
population but different parents was also 
successfully implemented for Meretrix meretrix 
culture (Lu ef a/ 2011). Due to the intensive 
mtDNA marker application, the universal 
primer of COI had been created based on the 
conserve gene sequence across the species to 
facilitate COI analyzing region on different 
metazoan species (Folmer ef a/. 1994). 

The study of Mollusca started from the study 
of Morris & Purchon (1981), their approach was 
to enlist and taxonomic study of available 
species; however, the story of mollusk study 
goes on, and Malaysian scientists and 
academician studied a different aspect of 
Mollusca; from taxonomy to diversity, along 
with the genetic aspect of this community 
(Hamli e¢ a/ 2012a, b; 2013; 2015; 2017; 2019; 
2020a, b, c; Foon ef af 2017; Al-Asif ef af 2020; 
AlL-Asif et al. 2021; Idris e¢ a/. 2017a, b). Genetic 
and molecular study of class Bivalvia and 
Gastropoda is not new in Malaysian territory, 
including East Malaysia (Sabah and Sarawak 
Provinces). Some previous study revealed the 
genetic diversity of Asian green mussel (Perna 
viridis) in the Waters of Sabah (Lau e¢ a/ 2018), 
DNA study of Neritid (Chee & Mohd Nor 
2016), genetic diversity of razor clam in Kuala 
Selangor, Malaysia (Hassan & Kanakaraju 2016). 
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There is no genetic approach of Venus clam 
(Meretrix spp.) was found in Malaysia, and the 
genetic relationship among the available Venus 
clam species is still unknown. Therefore, the 
present study aims to differentiate among 
Meretix spp. from different through the 
phylogenetic tree. 


MATERIALS AND METHODS 


Sample Collection 


A total of 30 samples of Meretrix hyrata, M. 
meretrix, and M. lusoria were collected from two 
areas, Kuching and Kabong, Sarawak Province, 
Malaysia (Fig. 1). Adductor muscle tissues from 
the females were separated from the shell and 
preserved with 95% of ethanol (Wang ef a/. 
2010) for further analysis at the Parasitology 
Laboratory, Department of Animal Science and 
Fishery, Universiti Putra Malaysia, at the Bintulu 
Sarawak Campus. 


DNA Extraction 


About 10 to 20 mg of adductor muscle tissue 
from each sample was used for extraction 
using a tissue DNA extraction kit (Vivantis) 
following the manufacturer's instructions and 
stored at -20 °C. The extracted DNA was then 
quantified through agarose gel electrophoresis. 


DNA Quantification and Gel Electrophoresis 


Extracted DNA was quantified following 
Rengarajan e¢ al. (2002) to ensure the amount of 
DNA within 1 wl. Each quantified DNA was 
adjusted to 100 wL to ensure the total amount of 
DNA sample used for amplification was equal. 
Quantification of genomic DNA was compared 
based on Lambda Hind III marker (23130 bp). 
Gel agarose concentration used for DNA 
quantification was 0.8% in 1X TAE (Trtis- 
Acetate-EDTA) buffer and the electrophoresis 
procedure was conducted following Lee ef a/. 
(2012). A total of 5 uL of extracted genomic 
sample and 1 ul of dye (EZ-Vision One dye) 
were added into the well. After all the samples 
and markers were loaded, electrophoresis was 
run for 40 minutes with 120 Volts of electricity. 
The running of electrophoresis was monitored 
until the red dye was moved 3/4 of the gel 
length, and the electrophoresis was stopped. 
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Figure 1 Study area of Meretrix spp. for genetic study 


The gel was then removed from the tank and 
was photographed under a UV light gel imager 
(Red™, Alpha Innotech). Quantification process 
based on appeared bands on the gel to be 
viewed with Alpha View software and band 
intensity from extracted samples was compared 
with Lambda Hind HI marker for 100 ng of 
DNA. 

While 1.5% of gel agarose concentration was 
used for the Polymerase Chain Reaction (PCR) 
through 0.3 g of agarose powder into 20 mL of 
1X TBE (Tris-Borate-EDTA) buffer, the 
mixture switled and heated in the 
microwave to ensute agarose powder was 
dissolved. Heated mixture was monitored to 
make sure the mixture was not boiled. After 
agarose powder dissolved in the buffer, the 
mixtute was left to cool and then was poured 
slowly into the prepared gel tank while ensuring 
that no bubble was trapped. Immediately, the 
comb was inserted into the poured mixture, 
while letting the agarose mixture to be cool and 
solidify for 30 minutes. After the gel became 
solid, the comb was removed and the running 
buffer 1X TBA or 1X TBE was poured into the 
gel tank until all the gel parts submerged. 


was 


Polymerase Chain Reaction (PCR) 
Amplification 


Amplification of DNA was performed in 
50 pL of PCR mixture using universal marker 
cytochrome C oxidase subunit I (COT) with 


sequence LCO1490: 5’--GGTCAACAAATCATA 
AAGATATTGG-3’? and =HCO2198: —5’- 
TAAACTTCAGGGTGACCAAAAAATCA-3’ 
(Folmer et ai 1994). Each PCR teaction 
consisted of 0.1 mM of reverse and forward 
primers, 0.2 mM of dNTP, 2.0 mM MeCh, 1x 
buffer, 1.25 units of Taq polymerase, 100 ng of 
template DNA, and ultrapure water to complete 
50 nL. 

A total of 50 wL of PCR mixture were 
amplified in XP Thermal Cycler Block (Bioer 
Technology co. LTD) with pre-denaturation at 
95 °C for 5 minutes, then followed by 35 cycles 
of 1 minute at 95 °C for denaturation, 1 minute 
for annealing at 40 °C, 1 minute 30 second at 
72 °C for extension and 72 °C for 7 minutes for 
final extension (Folmer e¢ ai 1994). The last 
cycle temperature was maintained at 4 °C. The 
PCR product then underwent electrophoresis. 


Gel Electrophoresis 


Gel agarose concentration used was 1.5% in 
1X TBE (Tris-Borate-EDTA) buffer. The 
electrophoresis procedure was carried out 
following Lee ef a/, (2012). A total of 100 ng of 
PCR product was added into the well with 1 pL 
of dye (EZ-Vision One dye) and 1 wh of 100 bp 
ladder (Promega) with 1 uL of dye (EZ-Vision 
One dye) at the end of the well. After all 
samples and markers were loaded, 
electrophoresis was run for 60 minutes with 80 
Volts of electricity. The running = of 
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electrophoresis was monitored until the red dye 
was moved 3/4 of the gel length, and the 
electrophoresis was stopped. The gel was then 
removed from the tank and was photographed 
under a UV light gel imager (Red™, Alpha 
Innotech). The appeared band was then com- 
pared to the ladder to identify the fragment size. 


Purification and Sequencing 


PCR products were purified before 
sequencing. According to the manufacturer's 
instruction, purification has been done to 
remove excess dinucleotide _ triphosphate 
(ANTP) using Themo Scientific GeneJet PCR 
Purification Kit. Afterward, purified samples 
wete sent to the First Base Laboratory, Malaysia, 
for DNA sequencing. Both strands of PCR 
products and primers were sequenced using the 
Applied Biosystems BigDye Terminator v3.1 
cycle sequencing kit. 


Statistical Analysis 


The obtained sequences were inspected and 
assembled using BioEdit version 7.0. The 
assembled sequence of M. /yrata, M. meretrix, and 
M. lusoria were then compared with all species of 
Meretrix spp. in GenBank of National Center for 
Biotechnology Information (NCBI) using Basic 
Local Alignment System Tool (BLAST) 2.2.31 
(National Center for Biotechnology Information, 
Bethesda, MD, USA [http://www.ncbi.nlm. 
nih.gov/BLAST/]). The highest — similarity 
sequence from the GenBank with Meretix spp. 
was selected for phylogenetic analysis in the 
present study. 


The present Meretrix spp. sequence, 16 
accessions sequence mtDNA, COI sequence, 
and 2 species with COI sequence as an outgroup 
which Mercenaria mercenaria and Paphia gallus 
(Table 1) in Genbank NCBI were aligned based 
on multiple sequences comparison by log- 
expectation using MUSCLE (Edgar 2004). 
Phylogenetic relationships were developed using 
Maximum Likelihood (ML) statistical methods 
with MEGA 6.0 statistical software (Tamura 
et al. 2013). 

The evolutionary history was inferred using 
the Maximum Likelihood method based on the 
Tamura 3-parameter model (Tamura 1992). The 
non-uniformity of evolutionary rates among 
sites was modeled using discrete Gamma 
distribution (+G) with 5 rate categories. The 
bootstrap agreement tree inferred from 500 


replicates (Felsenstein 1985) was taken to 
represent the evolutionary history of the 
analyzed taxa (Felsenstein 1985). Branches 


equivalent to partitions reproduced in more than 
50% bootstrap replicates were illustrated. The 
percentage of replicate trees in which the 
associated clustered together in the 
bootstrap test was shown next to the branches 
(Felsenstein 1985). The heuristic search's initial 
tree was obtained automatically by applying 
Neighbor-Join and BioNJ algorithms to a matrix 
of pairwise distances estimated using the 
Maximum Composite Likelihood (MCL) 
approach and then selecting the topology with a 
superior log-likelihood value. Evolutionary 
analyses were conducted in MEGAO6 (Tamura 
et al. 2013). 


taxa 


Table 1 List of Meretrix accession sequences from the GeneBank (NCBI) 


Meretrix accessions Geographical locations GenBank Acession Number Citation 

Meretrix lyrate Sarawak, Malaysia Present study Present study 

M. lusoria Sarawak, Malaysia Present study Present study 

M. meretrix Sarawak, Malaysia Present study Present study 

M. hrata China KC832317.1 Wu et al. (2014) 

M. hrata China J N898944.1 Cheng ef al. (2013) 

M. hrata China HM124580.1 Chen ef al. (2011b) 

M. lusoria China GQ903339.1 Wang et al. (2010) 

M. lusoria Japan AB853870.1 Yamakawa & Imai (2013) 
M. lusoria Japan AB853865.1 Yamakawa & Imai (2013) 
M. meretrix China GQ463598.1 He e¢ al. (2011) 

M. meretrix China HM124578.1 Chen ef af. (2011b) 

M. meretrix China JN043623.1 Wang et al, (2011) 

M. petechialis China EU145977.1 Ren et af. (2009) 

M. petechialis China HM124582.1 Chen e¢ af. (2011b) 
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Table 1 (Continued) 


M. petechialis China 
M. lamarckii China 
M. lamarckii China 
M. lamarckii Japan 

M. casta India 

Mercenaria mercenaria * Canada 
Paphia gallus * China 


HQ703177.1 Chen ef al, (2011a) 
GU071281.1 Wang ef al. (2011b) 
HM124579.1 Chen e¢ al. (2011b) 
AB059420.1 Hamaguchi ef a/. (2001) 
JQ773441.1 Ranjith e¢ ad. (2012) 
HM884239.1 Layton ef al. (2014) 
HQ703232.1 Chen ef al, (20112) 


Note: * = Outgroup. 


RESULTS AND DISCUSSION 


The COI marker produces a fragment size at 710 bp for M. Ayrata, M. meretrix, and M. Jusoria (Fig. 2). 
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Figure 2 Fragment size produced by Meretrix spp. using COI marker 
Notes: 1 = M. Arata, 2 = M. lusoria; 3 = M. meretrix. 


Basic Local Alignment System Tool (BLAST) 
analysis indicated that the M. /yrata samples from 
the present study had 99% similarity to the three 
accessions sequence (KC832317.1, JN898944.1, 
HM124580.1) of M. Arata compared to other 
accessions sequence (Table 2). Different 
accession sequences were only showed lower 
than 87% of similarities with M. Arata from the 
Sarawak locality. The present study on M. /usoria 
showed higher similarities with M. Arata with 
accessions sequence KC832317.1 and 
HM124580.1, and M. meretrix with accessions 
sequence HM124578.1 and JN043623.1. Other 
accessions sequence were less than 89% of 
similarities with the present study on M. /usoria. 
There was no similarity of more than 89% 
within M. meretrix from Sarawak Province with 


GeneBank NCBI accessions sequence of M. 
meretrix. "The higher similarity sequence was 
represented by the accession sequence 
HM124580.1 of M. Arata with 89% similarity. 
The nucleotide sequence for M. /yrata in the 
present study showed potential similarity with 
the accession nucleotide sequence of M. Arata 
from the GenBank NCBI with more than 90% 
similarity compared to M. meretrix and M. /usoria. 
Some previous studies conferred the 
morphology feature on M. Arata as being the 
easiest to distinguish from M. meretrix and M. 
lusoria (Hamli e¢ al 2012b, 2015, 2016, 2017). 
However, the nucleotide sequence for M. 
meretrix and M. lusoria showed lower than 90% 
similarity with the accession nucleotide sequence 
from the GenBank NCBI. Hence, the BLAST 
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analysis could not confirm and distinguish M. 
meretrix and M. /usoria from the Sarawak locality 
through the COI sequence. Current findings 
wete contradictory with results of the 
morphology and morphometric studies, which 
distinguished differences between M. meretrix 
and M. /usoria. However, Wang & Dunbrack 
(2004) reported that the application of BLAST 
only produces short and conserve fragment 
alignment. Thus, the analysis was frequently 
unable to identify many relationships, 
particularly below 40% sequence identity. Low 
sequence identity of fragment creates 
uncertainty of gaps and insertion, thus 
compromising confidence in overall modeling of 
protein related to genome sequence annotation 
(Gan et al. 2002). M. lusoria and M. meretrix were 
hard to identify by rough morphological 
inspection due to similar features, particularly 
the outer shell. However, a detailed examination 
of the inner shell generated morphological 
variation of the distinct profile of the pallial 
sinus scat between M. /usoria and M. meretrix. 
Yet, using the BLAST application to confirm 
the differences between the two species 
appeared to provide uncertain results. 

Rather than depending on sequence identity 
percentage, phylogenetic relationship perhaps 
provides mote justified evidence on the present 
genome characteristic of Meretrix spp. to be 
evaluated against the accession sequence in the 
GenBank. The current phylogenetic tree was 
constructed based on Maximum Likelihood 
(ML), which applied an algorithm that calculates 
the probabilistic approaches to build a 


relationship tree between compared species 
based upon the nucleotides or amino acid 
sequences (Guindon & Gascuel 2003). Meretrix 
spp. from the Sarawak locality indicates the 
close relationship between M. /usoria and M. 
meretrix. The variation in the COI sequence of 
this Meretrix spp. can be distinguished, hence, we 
are able to clarify the intricacy conferred using 
the BLAST application. 

Moreover, M. /usoria and M. meretrix had 
similar ancestors represented by M. yrata from a 
similar locality. Therefore, COI sequence in M. 
yrata had a low level of sequence identity among 
Meretrix spp. This was in agreement with Chen e¢ 
al. (2009) and Wu e¢ a/. (2014), which suggested 
that M. Arata has the level of sequence that 
portrays as the ancestor to other Meretrix spp. by 
referring to COI and transfer RNA (tRNA), 
respectively. 


Sequence Characteristics 


A total of 21 nucleotide sequences were 
analyzed from three Meretrix spp. of Sarawak 
locality, 16 accessions sequence of Meretrix spp. 
and two outgroup species represented by 
Mercenaria mercenaria and Paphia gallus. Codon 
positions included were 1% + 2% + 3% + 
Noncoding. All positions containing gaps and 
missing data were eliminated. There were a total 
of 435 positions in the final dataset. The 
nucleotide base composition sequences were 
comprised with mean of Thiamine (TI) 37.9%, 
Cytosine (C) 15.4%, Adenine (A) 27.4%, and 
Guanine (G) 19.4% (Table 3). 


Table 2 Similarity percentage of Meretrix spp. in the present study with 16 selected Meretrix spp. mtDNA accessions 


sequence from GenBank (NCBI) 


Accession sequence Speci M. hyrata M. lusoria M. meretrix 
GenBank PROee (%) (%) (%) 
KC832317.1 M. hrate 99 90 88 
JN898944.1 M. byrata 99 < 83 < 84 
HM124580.1 M. hrata 99 90 89 
GQ463598.1 M. meretrix 86 86 85 
HM124578.1 M. meretrix 87 90 87 
JN043623.1 M. meretrix 87 90 87 
GQ903339.1 M. lusoria < 85 86 86 
AB853865.1 M. lusoria 85 87 85 
AB853870.1 M. lusoria < 85 86 85 
EU145977.1 M. petechialis 86 85 85 
HM124582.1 M. petechialis < 85 85 85 
HQ703177.1 M. petechialis < 85 86 85 
GU071281.1 M. lamarckii < 85 < 83 86 
AB059420.1 M. lamarckii < 85 < 83 87 
HM124579.1 M. lamarckii < 85 < 83 86 
JQ773441.1 M. casta < 85 88 87 
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Table 3 Nucleotide base composition for accession sequences 


Species Accession T% C% A% G% Total 

Meretrix lyrata Sarawak 24.8 19.5 42.1 13.6 435.0 

M. lusoria Sarawak 24.4 20.2. 41.8 13.6 435.0 

M. meretrix Sarawak 23.7 20.7 42.1 13.6 435.0 

M. byrata KC832317.1 41.8 14.0 23.7 20.5 435.0 

M. byrata JN898944.1 42.1 13.8 24.1 20.0 435.0 

M. byrata HM124580.1 42.1 13.8 24.1 20.0 435.0 

M. lusoria GQ903339.1 42.1 13.6 23.0 21.4 435.0 

M. lusoria AB853870.1 23.0 20.5 43.0 13.6 435.0 

M. lusoria AB853865.1 22.3 22.1 43.0 12.6 435.0 

M. meretrix GQ463598.1 42.8 13.1 23.0 21.1 435.0 

M. meretrix HM124578.1 41.6 14.3 22.3 21.8 435.0 

M. meretrix JN043623.1 42.1 13.8 24.1 20.0 435.0 

M. petechialis EU145977.1 42.8 13:1 23.0 21.1 435.0 

M. petechialis HM124582.1 41.8 13.8 23.7 20.7 435.0 

M. petechialis HQ703177.1 42.3 13.3 23.7 20.7 435.0 

M. lamarckii GU071281.1 43.7 13.8 19.8 22.8 435.0 

M. lamarckii HM124579.1 43.7 13.6 20.7 22.1 435.0 

M. lamarckii AB059420.1 42.3 13.8 22.1 21.8 435.0 

M. casta JQ773441.1 42.3 14.7 23.2 19.8 435.0 

Mercenaria mercenaria HM884239.1 40.8 14.7 21.1 23.4 436.0 

Paphia gallus HQ703232.1 44.4 12.0 23.0 20.7 435.0 

Average 37.9 15.4 27.4 19.4 435.0 
Phylogenetic Analysis the China locality were clustered with M. meretrix 


A discrete Gamma distribution was used to 
model evolutionary rate differences among sites 
(5 categories (+G, parameter = 2.0104)). 
Maximum Likelihood (ML) analysis conducted 
using the Tamura 3-parameter model has 
constructed five major clades on Merefrix spp. in 
the present study with designed out-groups 
Mercenaria mercenaria and Paphia gallus. These 
major clades were closely related to each other 
at the 50% evidence of bootstrap, which 
grouped as Meretrix genus (Fig. 3). Two 
monophyly groups were formed, represented by 
Clade 3 and 4, while Clade 1, 2, and 5 were 
formed, paraphyly groups. Meretrix spp. from 
the Sarawak locality was clustered in Clade 1 
alongside two M. /usoria of Japan locality with 
accession sequence number AB853865.1 and 
AB853870.1. Clade 2 comprised two M. meretrix 
from the China locality with accession sequence 
number HM124578.1 and JN898944.1. Clade 3 
consisted of M. /yrata from the China locality 
with accession sequence number KC832317.1, 
HM124580.1, JN898944.1, and M. casta from the 
India locality with accession sequence 
JQ773441.1. Meretrix lamarcku from the Japan 
locality (AB059420.1) was clustered with M. 
lamarckii from the China locality (GU071281.1; 
HM124579.1) in Clade 4. All M. petechialis 
(HM124582.1, HQ703177.1, EU145977.1) from 


(GQ463598.1) and M. /usoria (EU145977.1) 
from the China locality, respectively, in Clade 5. 

Phylogenetic analysis on three Meretnix from 
the Sarawak locality is closely related to the 
M. lusoria from the Japan locality. However, the 
result contrasts with M. Arata from BLAST 
analysis which inferred a high similarity 
sequence with M. /yrata accession sequence. The 
application of ML analysis for the phylogenetic 
tree was calculated using the substitution rate in 
the sample sequence and accession sequence, 
which produced different inferences than 
BLAST analysis (Guindon & Gascuel 2003). 
Furthermore, Meretrix spp. from the Sarawak 
locality and M. /usoria from the Japan locality had 
a high percentage of Cytosine and Adenine than 
that of other Meretrix spp. from different 
localities. Meretrix spp. from the Sarawak locality 
only had a relationship with other Meretx spp. 
localities at more than 50% of bootstrap scoring, 
which is inferred as under Meretrix genus. 
Inconsistent cladistic group between M. meretrix 
(GQ463598.1) M. lusoria (GQ903339.1), M. casta 
(JQ773441.1), and M. petechialis (HM124582.1, 
HQ703177.1, EU145977.1) from the present 
study indicated that the three species had a close 
relationship following Wu ef a/. (2014) and Chen 
et al. (2009). Hence, Chen e¢ a/. (2009) suggested 
that M. /usoria and M. petechialis as a junior 
synonym to M. meretrix. 
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Meretrix petechialis HM124582.1 


a Meretrix petechialis. — HQ703177.1 
5 9 ‘Meretrix lusoria GQ903339.1 
‘Meretrix petechialis EU145977.1 

‘Meretrix meretrix GQ463598.1 

‘Meretrix lamarckii AB059420.1 
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Paphia gallus HQ703232.1 


Figure 3 Phylogenetic tree of Merefrix accession sequence inferred from the Maximum Likelihood analysis using 


Tamura three-parameter model 


Note: Only bootstrap score greater than 50% are shown. 


The present study on Meretrix spp. from the 
Sarawak locality was able to differentiate COI 
sequences among M. Arata, M. meretrix, and 
M. lusoria. M. lusoria was closely related to 
M. meretrix with strong bootstrap supporting 
evidence at 96% scoring. M. Arata became the 
ancestor to M. /usoria and M. meretrix with strong 
bootstrap evidence at 100% of the scoring. 
Meretrix spp. from the Sarawak locality is closely 
related to the M. /usoria of the Japan locality with 
100% scoring on bootstrap supporting evidence 
to form Clade 1. Meretrix spp. from the Sarawak 
locality was closely associated with other 
accession sequences for Meretrix spp. at 
bootstrap supporting evidence at 50% scoring. 

Based on ML phylogenetic analysis, all 
M. Wrata with accession sequence KC832317.1, 
JN898944.1, and HM124580.1 formed the 
cladistic group with strong bootstrap evidence at 
99%. All M. lamarckii with accession sequence 
GU071281.1, AB059420.1, and HM124579.1 
also formed cladistic at 94% of bootstrap 
evidence. M. casta (JQ773441.1) formed an 
individual clade which indicated this species 
distance from other species except for M. /yrata. 
Other Meretrix spp. accession sequence, such as 
MM. meretrix (GQ463598.1) M. = lusoria 
(GQ903339.1), and M. petechialis (HM124582.1, 
HQ703177.1, EU145977.1) formed inconsistent 
cladistic groups, while MM. meretrix with accession 
sequence HM124578.1 and JN043623.1 formed 
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individual clade apart from other M. meretrix 
(GQ463598.1) at 99% evidence of bootstrap 
scoring. 

Mitochondria DNA is 
higher transition mutation rates, particularly 


known to have 


synonymous sites in the nucleotide sequence 
that affected the rates of evolution in 
organisms (Brown ¢ al 1982; Caterino ef al. 
2000; Overton & Rhoads 2004). This fast rate of 
evolution is one of the features suitable for 
phylogenetic study, particularly in higher taxa. 
Moreover, higher taxa less affected by the 
recombination (Guo ef a/ 2006), uniparental 
(Passamonti e¢ al 2003; 
Passamonti & Plazzi 2020), heteroplasmy with 
paternal leakage (Bromham ¢éf a/, 2003; Wolff e¢ 
al, 2013; Mastrantonio et a/ 2019), and selective 
sweeps (Ballard & Rand 2005; James e¢ a/, 2016; 
Hill 2019). Despite proper application on higher 
this is 
sequences at the species level take into account 
those which ate more defected. The evolution 
rate is thermal adaptation, 
mitochondrial, nuclear interaction, and infection 
with Wolbachia spp. (Ballard & Rand 2005). 
Hence, Meretrix spp. from the Sarawak locality 
possibly has a different evolution rate in 
mtDNA, which caused divergence from the 
Meretix spp. from other localities, particularly 
from the China regions. 


inheritance 


taxa, the contrast when nucleotide 


influenced by 
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The current phylogenetic analysis is difficult 
to verify that species Meretix spp. belong to the 
Sarawak locality. Bivalves had high variability in 
size and gene group of mtDNA, which vary 
extensively between species from a similar genus 
(Xu ef al. 2012). Meretrix spp. from the Sarawak 
locality are separated in terms of geographical 
distance. The geographical distance affects the 
mtDNA feature as described by Mishmar ef a/. 
(2003) that 
restricted by the geographical distance with no 
mixing of mtDNA haplogroups from the 
different geological region. Therefore, other 


maternal inheritance can be 


environmental conditions in each region 
influence the organisms' mutation rate 
(Massey & Buckling 2002). Ambiguous 


Meretix spp. identification from the Sarawak 
locality might be possibly related to other several 
features, such as sample used. 

The mtDNA marker feature is the maternal 
which supplies the 
information in nuclear markers (Okumuu & 
Ciftci 2003). Therefore, the female individual's 
tissue selection is crucial during the genetic 


inheritance missing 


study to ensure reliable information regarding 
the species ancestor. However, female selections 
for the genetic 
mentioned in each previous study for the 
GenBank the 
doubtness reliability of accession 


reseatch have not been 


accession 
of the 
sequence in the GenBank, leading to the 


sequence, causing 


divergence of Meretrix spp. from the Sarawak 
species 
sequence. Moreover, the accession sequence is 
doubted due to the 
morphological identification for each species 


locality compared to the previous 


well non-intensive 
which only depending on the rough observation. 
The present study suggested that appliance on a 
single marker is insufficient to support species 
identification and phylogenetic 
analysis. Therefore, additional variance types of 


species-level 


molecular markers such as allozyme and nuclear 


matkers are required to create substantial 
evidence on species verification to ensure the 
study's —_ genetic 


Depending on the genetic 


intensive 
approachability. 
approach, the lack of intensive morphological 


morphological 


evidence would create ambiguous species 


naming in taxonomy. 


CONCLUSION 


The molecular approach application was able 
to distinguish three Meretrix spp. from Sarawak 
Province, Malaysia. Moreover, M. Jyrata is 
inferred to as the ancestor of M. meretrix and 
M. lusoria from Sarawak. Regardless of ancestor 
identification and species distinctions, the COI 
gene in the present study was unable to verify 
which species Meretrix spp. belonged to the 
Sarawak locality. Several disadvantages to the 
COI sequence that influence the phylogenetic 
analysis involve the consistent accession 
sequence from GenBank as the main reference 
for the current study. The present study 
suggested that additional DNA markers for 
species identification and phylogenetic analysis 
are recommended in supporting the intensive 
morphological identification. 
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